% GAUTIER LE BIHAN - 2020
% Replication files for "Shocks vs Menu Costs: Patterns of Price Rigidity in an Estimated Multi-Sector Menu-Cost Model?" Review of Economics and Statistics
%
% This code produces  (part of)  Table D
clear;
tic
addpath('..\..\Utilities')  

load actual_moments_k

tic
cd ..\..\Simulations_VC\Counterfactual_sim\MS_produits_KR_MC

load ..\..\Estimation_param\MS_produits_KR_MC\actual_moments_k
weight=actual_moments_k(:,12);
sum(weight)
secteur=actual_moments_k(:,1);

load ..\..\Simulations_VC\MS_produits_KR_MC\stat4
param=stat4;


meanw_data_all=weight'*param(:,5:9)/sum(weight);
for i=1:5
stdw_data(i)=sqrt(weight'*((param(:,i+4)-meanw_data_all(i)).^2)/sum(weight));
end

stdnw_data=std(param(:,5:9));
meannw_data_all=mean(param(:,5:9));

[stdnw_data./meannw_data_all; stdw_data./meanw_data_all]


%HE
meanw_data_he=weight(secteur ~=5)'*param((secteur ~=5),5:9)/sum(weight(secteur ~=5));
for i=1:5
stdw_data_he(i)=sqrt(weight(secteur ~=5)'*((param((secteur ~=5),i+4)-meanw_data_he(i)).^2)/sum(weight(secteur ~=5)));
end
stdnw_datahe=std(param((secteur ~=5),5:9));
meannw_data_he=mean(param((secteur ~=5),5:9));

[stdnw_datahe./meannw_data_he ; stdw_data_he./meanw_data_he]





load ..\..\Simulations_VC\MS_produits_KR_MC\stat_simu_out.mat

for jj=2:4
    
    param=stat_simu_out(jj).stat_simu_out(:,6:10);
    
    meanw_data=weight'*param/sum(weight);
for i=1:5
std_w(i)=sqrt(weight'*((param(:,i)-meanw_data(i)).^2)/sum(weight));
end

stdw(jj,:)=std_w;
mw(jj,:)=meanw_data;

stdnw(jj,:)=std(param);
mnw(jj,:)=mean(param);
%HE
mean_data=weight(secteur ~=5)'*param((secteur ~=5),:)/sum(weight(secteur ~=5));
for i=1:5
std_w(i)=sqrt(weight(secteur ~=5)'*((param((secteur ~=5),i)-mean_data(i)).^2)/sum(weight(secteur ~=5)));
end

stdw_he(jj,:)=std_w;
mw_he(jj,:)=mean_data;

stdnw_he(jj,:)=std(param((secteur ~=5),:));
  mnw_he(jj,:)=mean(param((secteur ~=5),:));
end

fprintf('Weighted');

fprintf('ALL');
[stdw_data./meanw_data_all ; stdw./mw]'
fprintf('non energy');
[stdw_data_he./meanw_data_he ; stdw_he./mw_he]'

fprintf('UnWeighted');

fprintf('ALL');
[stdnw_data ; stdnw]'
fprintf('non energy');
[stdnw_datahe ; stdnw_he]'



